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ABSTRACT 



We report on all the INTEGRAL and Swift data collected during the first outburst observed from IGR J18179-1621 . The broad-band 
spectral analysis showed that the X-ray emission from the source is heavily absorbed (_/V H -10 23 cm -2 ), and well-described by a flat 
power-law with a high energy rollover (cutoff energy 9-12 keV, e-folding energy 4-7 keV). We found some evidence of a cyclotron 
absorption feature at 22 ± 1 keV. Together with the pulsations at 1 1.8 s discovered in the XRT data, this evidence would suggest that 
IGR J18179-1621 is an obscured, magnetized, accreting neutron star that is possibly part of a supergiant high-mass X-ray binary or a 
Be X-ray binary system. 
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1. Introduction 

IGR J18179-1621 was discovered by the hard X -ray imager 
ISGR I on-board INTEGRAL on 2012 February 29 dTiirler et al.1 
|2012|) . The estimated source flux in the 20-40 keV energy band 
was 16+1 mCrab, (1.23 + 0.08) x lO" 10 erg cm -2 s _1 . Follow- 
up observations with Swift and Chandra permitted to obtain 
the refined position of the source at aj2000=18 h 17 m 52 s .19 and 
5j2ooo=-16°21'31'.'7 (error-circle radius of 07.6 at 90% confi- 
dence level), and led to the identification of the cand idate in- 
frared (IR) counterpart 2MASSJ18 175218- 16213 16 dLi et alj 
I20TI iPaizis eTaIll2012h . The narrow-field instrument on-board 
Swift, XRT, also measured a high absorption column-density in 
the direction of the source (~10 23 cirT 2 ) and det ected pulsa- 
tions in its X-ray em ission at a period of 11.82 s dTiirler et alj 
l2012_lHalpernl2012l). whi ch was confirmed using FERMI/GBM 
(Finger & Wils on-Hodgell2012l) . This evidence suggested that 
IGR J18 179-1621 is a highly obscured X-ray pulsar, possibly 
part of a high-mass X-ray binary system. 

In this paper, we report on the analysis of all the 
INTEGRALflSGRl and Swift /XRT target-of-opportunity ob- 
servations we requested to be performed in the direction of 
IGR J18179-1621, and discuss different interpretations of the na- 
ture of the source. 



2. Data analysis and results 

2.1. INTEGRAL 

IGR J 1 8 1 79- 1621 was detected simultaneously by IB IS/ISGRI 
(15-500 keV. IUbertini et al.ll2003ULebmnet al. 2003) and JEM- 
X (3-25 keV. iLund et al.ll2003l) for the first time during satel- 
lite revolution 1145, and observed throughout revolutions 1146 
and 1147 (from MJD 55 986.1 to MJD 55 992.5). Further ob- 
servations were inhibited by a bright solar flare that forced the 
instrument on-board INTEGRAL to enter in "safe-mode" with 



Table 1. INTEGRAL observation log and effective exposure 
times of the spectra of IGR J 1 8 1 79- 1 62 1 . 



Rev. 


START TIME 


STOP TIME 


ISGRI 


JEM-X1 


JEM-X2 




(MJD) 


(MJD) 


(ks) 


(ks) 


(ks) 


1145 


55986.1 


55987.8 


55.6 


23.7 


23.1 


1146 


55988.9 


55990.8 


62.8 


35.3 


34.6 


1147 


55991.8 


55992.5 


37.4 


21.8 


20.9 


Note 


A selection for 


source off-set angles 


has been applied, see 


Sect.im 



no scientific data available^ A complete log of these observa- 
tions is provided in Table Q] All INTEGRAL data were ana- 
lyzed using version 9.0 of the OSA software distributed by the 
ISDC dCourvoisier et alj 120031) . INTEGRAL observations are 
commonly divided into "science windows" (SCWs), i.e. point- 
ings with typical durations of ~2-3 ks. Because of uncertainties 
in the instrument responses for high off-axis angles, we selected 
for the spectral analysis only those SCWs for which the off-set 
angle of the source was <3.5 deg and <12 deg, respectively, for 
JEM-X and ISGRI. The total effective exposure-time resulting 
after this selection is reported in Table [1] In all revolutions, the 
source was relatively faint for the JEM-X sensitivity, thus we 
performed a spectral extraction by using the option flag—1 in 
the local source catalog. We checked a posteriori that this did 
not affect the discussion on the main spectral results reported in 
Sect. 12.21 Owing to the relatively low signal-to-noise ratio (S/N) 
of the INTEGRAL data, no meaningful timing analysis could be 
carried out. 



Send offprint requests to: enrico.bozzo@unige.ch 



1 See http://www.isdc.unige.ch/integral/operations/reports 

2 This means that the source spectrum is extracted 
at the given catalog position even if the source is 
not detected within a single science window, see 
http://www.isdc.unige.ch/integral/download/osa/doc/current/osa_um 
_jemx/node40.html 
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Table 2. Swift observations log of IGR J18 179-1621. 



OBS" 


START* 


STOP 4 Exp N u ' 


r 


fob/ xljd.o.f. 




(day) 


(day) (ks) 




(C-stat/d.o.f.) 



Note: all data were fit with a simple absorbed power-law model. 

a : Only the last two digits of the observation ID. 000322930** are shown. 

b : Start and stop times of the observations are from MJD 55 986, 29 Feb 2012. 

c : The absorption column density is in units of 10 22 cm" 2 . 

d : The absorbed X-ray flux is in units of 10~ 10 erg cm -2 s -1 

e : Owing to the relatively low statistics, the three observations 15, 16, and 17 were 
merged to allow for a meaningful spectral extraction. 

IGRJ1 8 179-1621 (Swift/XRT) 



"1145 1146 1147 
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Fig.l. Swift /XRT lightcurve of IGR J18179-1621 during the 
~22 days monitoring reported in this paper (time bin is 1 ks). The 
time intervals covered quasi-simultaneously by INTEGRAL is 
also reported with a solid line. 



2.2. Swift 

Swift /XRT observations were performed from a few hours up to 
~22 days following the discovery of the source (29 Feb 2012- 
22 Mar 2012, see Tabled. Swift/XRT data were collected only 
in photon-c ounting (PC) mode a nd analyzed by using standard 
procedures (Burrows et al. 2005). The XRT data were processed 
with the xrtpipelin e (v.0. 12.6) and the latest calibration files 
available (see also iBozzo et al. 1 120091) . When required, we cor- 
rected PC data for pile-up, and used the xrtlccorr to account for 
this correction in the background-subtracted light curves. Source 
and background-event lists (time resolution 2.5 s) were barycen- 
tered by using the barycorr tool. 

All XRT spectra were well-described by a power-law model 
with a relatively large absorption, jVh^IO 23 cirT 2 (we used 
the phabs model in Xspec12), and flat photon index (r~0; see 
Table |2J. In order to obtain a better estimate of the source X- 
ray spectral properties, we performed a quasi-simultaneous fit 



of the XRT, JEM-X and ISGRI data. XRT observations from 
00032293001 to 00032293006 were summed-up together to per- 
form the quasi-simultaneous fit with the JEM-X and ISGRI data 
in revolutions 1 145 and 1 146, where the source displayed only 
a moderate change in the X-ray flux (we verified that no signif- 
icant spectral variation occurred in the two revolutions by per- 
forming separate fits to the data). As the source flux during rev- 
olution 1 147 is significantly lower than during the previous one, 
we performed a separate fit by combining the XRT spectrum ob- 
tained from the observation 00032293007 with the JEM-X and 
ISGRI spectra extracted during revolution 1 147. The results of 
this analysis are shown in Fig.|5]and reported in Table[3] 

Spectra from the first data-set (1 145+1 146) were fit first us- 
ing a number of phenomenological models, including in all cases 
a normalization constant to take into account the systematic un- 
certainty in the inter-calibrations among the instruments and the 
possible variability of the source (the constant was fixed to 1 
for XRT). A single absorbed power-law model failed to pro- 
vide an acceptable fit to the spectra Cfj ed -3.5). The addition 
of a cut-off at the higher energies (highecut*pow in Xspec) sig- 
nificantly improved the fit (^ 2 ed ^1.5), but left evident residuals 
above ~10 keV (see Fig. [5]). T he widely used NP EX and Fermi- 
Dirac cut-off models (see e.g. JCoburn et alJl2002h did not give a 
better fit (y 2 ed ^ 1 .7) and left very similar residuals. We also tried 
the most widely used Comptonization models: BMC comprises 
black-body (BB) radiatio n up-scattered by a spheric ally sym- 
metric in-falling plasma ( Shrader & Titarchuk 1998) and pro- 
vided a slightly better fit (x^ ed - 1 -4). However, the temperature of 
the thermal seed photons was significantly higher (~3 keV) than 
expecte d in the case of an accreting neutron star (0. 8-1.5 keV; 
see e .g., [Orlandini 2006). Similarly, a comptt model dTitarchukl 
1 19941) . which accounts for the thermal Comptonization of BB 
seed photons, was bounded to a temperature of the BB ^2.5 keV. 
In all models, evident residuals were left above ~ 1 k eV. As sug- 
gested by the preliminary data analysis reported bv lTiirler et all 
the combination of a cut-off power-law model plus a cy- 
clotron absorption line at ~20 keV gave a reasonably good fit to 
the data. However, the cyclotron feature parameters could only 
be poorly constrained; we thus fixed the line width (cr cyc i) to 
3 keV, which is a rea sonable value for an X-ray pulsar (see, e.g., 
iFerrigno etaT]|2009l and Table [3). In revolution 1147, the S/N 
of the data was relatively low compared to data in revolutions 
1 145 and 1 146, and the possible cyclotron feature was not for- 
mally required in the fit (an acceptable fit could be obtained by 
using an absorbed power-law model with a high energy cut-off, 
see Table [3). A further discussion of the interpretation of these 
results is given in Sect. [4] 

We performed a detailed timing analysis of the XRT data. 
Pulsations were detected at a period of ^11 .82 s in all obser- 
vations of a sufficiently long duration. In each of these ob- 
servations, we measu red the spin period by using Z 2 statistics 
dBuccheri et al.l[l983l) and determined the associated uncertainty 
using a Monte-Carlo approach, where we simulated 1000 event 
files with the average modulation of the data and repeated the 
Z 2 search on each of them. As detailed in the appendix, we es- 
timate the reliability of the observation by setting a threshold 
confidence level (c.l.) at 99%. Only four observations were char- 
acterized by a sufficient high quality to permit a reliable deter- 
mination of the spin period and its associated uncertainty. In the 
others, the spin period was determined at only a low confidence 
level and we were unable to reliably estimate the correspond- 
ing uncertainty (the folded lightcurve is reported for each obser- 
vation in Fig. [2] to allow an easier comparison). For each XRT 
observation, we also computed the fractional root mean square 
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Table 3. Results of the combined XRT/JEM-X/ISGRI spectra of IGR J18 179-1621. using an absorbed power-law with high-energy 
exponential roll-over and an absorption feature. 



Rev. 




r 


keV 


£fold 

keV 


Ecycl 

keV 


& cycl r cycl 

keV 


QsGRI 


Qemxi 


QEMX2 


^1-10 keV 

(10- 


*~4-20 keV 

erg cm"" 


^20-50 keV 




1145-6 


14-0!!;? 


0.5±0.2 


11.0±0.8 


6 2 +0 - 9 

u ' -0.6 


21.9!?? 


3.0* 8.1+1' 


1 2 +0 - 5 

1 ■•'-0.4 


1.0±0.1 


1.3±0.1 


2.6 


6.5 


1.2 


1.0/173 


1145-6" 


14-6!!:? 


0.5±0.2 


10.8±0.6 


t: 7 +0.6 
J - '-0.5 






1.0±0.2 


1.0±0.1 


1.3±0.1 


2.6 


6.6 


1.2 


1.5/175 


1147" 


n.i!g 


0.1±0.6 




4 8 +1 - 4 






0.8!^ 


1 4 +0 - 4 

i- -0.3 


1 4 +04 

1 -^-0.3 


1.6 


5.5 


0.6 


1.1/40 



a : The is in units of 10"- cm . 

. This parameter was fixed (see text for details). 
N : In these cases, no cyclotron feature was included in the fit. L Absorbed flux. 



Table 4. Spin period measurements from the XRT data. 
Uncertainties are at the Icr confidence level 
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OBS" 


START TIME 


DURATION 


P spi „ 


Hag" 


Fractional 




(MJD) 


(days) 


(sec) 




nns 


01 


55986.877 


0.037 


11.8230+0.0018 


Y 


0.20+0.03 


02 


55987.653 


0.012 


11.819+0.006 


Y 


0.21+0.02 


04' 


55988.957 


0.040 


11.7988+0.0016 


Y 


0.19+0.02 


05 


55989.395 


0.074 


11.7860 


N 


0.17+0.03 


06 


55990.632 


0.110 


11.8053 


N 


0.18+0.03 


07 


55992.261 


0.139 


11.8151 


N 


0.24+0.03 


08' 


55993.572 


0.247 


11.8271+0.0003 


Y 


0.29+0.03 


10 


55997.219 


0.008 


11.81 


N 


0.25+0.04 


12 


55999.367 


0.335 


11.8447 


N 


0.31+0.05 


15+16+17 


56006.550 


1.243 


11.81 


N 


0.46+0.07 



Only the last two digits of the observation ID.000322930** are shown. 
b : The flag "Y" indicates a spin-period determination that is reliable 
at more than 99% c.l., "N" otherwise. 
' : To be considered with caution, see AppendixlAl 



(rms) and its uncertainty using the same Monte Carlo approach 
described above (we simulated 1000 pulse profiles for each ob- 
servation, see Table No significant variation in the rms with 
luminosity was found. All measurements are compatible with an 
rms value of 20%. 

All the results of the timing analysis are summarized in 
Table 2] where we report the uncertainties associated with each 
spin-period measurement (l<x c.l.) and a flag expressing the 
reliability of the measurement. Reliable periods from XRT, 
together with previous determinations from Fermi/GBM (see 
Sect. [TJ, are plotted in Fig. [3] The first two determinations 
(Observations 00032293001 and 00032293002) gave a value of 
the spin period compatible with that determined from the GBM 
data, whereas the spin-period values measured from observa- 
tions 00032293004 and 00032293008 differed significantly from 
that. As discussed in Appendix lAl the origin of this deviation 
cannot be firmly established at present, thus these two measure- 
ments should probably be taken with caution. 

The XRT monitoring showed that the source began a rela- 
tively rapid decay in its X-ray flux after the initial brightening, 
and became fainter by a factor of ~40 in about 22 days. On 
MJD 56007.78, the source was too faint to obtain meaningful 
spectral measurements with the available exposure times of the 
XRT pointings (1-2 ks), thus no further observations were per- 
formed. 



3. Counterparts 

The arcsecond localization of IGR J18 179-1621 obtained first 
through the Swift observations and then refined with Chandra 
(see Sect. Q]), permitted to search for possible associated coun- 
terparts to the source i n the optical and IR domains. A search 
in the 2MASS catalog dSkrutskie et alJl2006l) revealed that the 
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Fig. 2. Normalized pulse profiles of the XRT observations (not 
phase-connected). The XRT observation number (ID) are indi- 
cated by the two last digits for each pulse profile. Each pulse has 
been divided by its average and vertically displaced by a con- 
stant value for clarity. Alignment in the figure is such that in 
each case the phase of the first Fourier component is 0.5. The 
solid lines are obtained by the Fourier decomposition using the 
first two components. 



object 2MASS 18175218 - 1621316 is the most likely IR coun- 
terpar t of the X-ray pulsar (see also iLi et al.ll2.012l : iPaizis et all 
120121) . We verified that the 2MASS object has been detected in 
K band with a magnitude of K = 13.14 + 0.04, while in the J and 
H bands only upper limits were reported (magnitudes of 16.7 
and 15.9, respectively). No obvious counterpart could be found 
at the position of IG R J 1 8 1 79- 1 62 1 from the USNO B-l catalog 
dMonet et al.l2003l) . At present, the lack of any detailed informa- 
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Fig. 3. Spin periods determined using the Z 2 statistics. Filled 
symbols are the subset of the determinations reliable at 99% 
c.l. obtained in the present work, open symb ols are the de- 
terminations of iFinger & Wilsor>H odge (2012) obtained using 
Fermi/GBM. 




Fig. 4. Images of the field of view around IGR J 1 8 1 79- 1 62 1 from 
the 2MASS archive (4'x4' ) Left: image in the J-band together 
with the Swift (dashed circle) and Chandra (solid circle) local- 
izations of IGR J18179-1621. The solid lines represent the lo- 
calization of the proposed IR counterpart to IGR J18179-1621, 
2MASS 718175218 - 1621316, which lies within the X-ray er- 
ror circles. Right: Same as previous figure but for the K-band, in 
which the counterpart is detected. 

tion in the J and H bands, combined with the unknown distance 
to the source, does not allow us to significantly constrain the na- 
ture of the companion to IGR J18179-1621. 

4. Discussion and conclusions 

IGRJ18179-1621 was detected for the first time with 
INTEGRAL during the period ranging from 55 986.1 MJD to 
55 992.5 MJD, even though Fermi/GBM data indicated that the 
outburst of the source might have already started on 55977 MJD. 
A follow-up monitoring with XRT was initiated shortly af- 
ter the discovery and continued until 56007.8 MJD, measur- 
ing a decay in the source X-ray flux of a factor of -40 (from 
2.9xl0~ 10 erg cirT 2 s"'to 8xl0~ 12 erg cirT 2 s _1 ) in -22 days. 
The detection of pulsations at ^1 1.82 s with the GBM and XRT, 
led soo n to classifyi ng IGR J18 179-1621 as an accreting X-ray 
pulsar (Halpern 2012|). 

The spectral analysis of all the available INTEGRAL and 
Swift data collected during the outburst revealed that the X-ray 
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Fig. 5. Our XRT/JEM-X/ISGRI combined spectra during the 
INTEGRAL revolutions 1 145 and 1 146 (see text for details). 
The best-fit model is obtained by using a power law with a high- 
energy cut-off plus a cyclotron absorption-line at ~20 keV. The 
mid panel shows the residuals from the fit when the cyclotron 
feature is not included. 



emission from the source was heavily obscured. We estimated 
an absorption column density of about an order of magnitude 
higher than th e Galactic value expected i n that direction (Nu = 
1.6xl0 22 cm~~ 2 lDickev & Lockma nl 19901) . In the brightest phase 
of the outburst (fromMJD 55 986.1 to MJD 55 990.8), the broad- 
band XRT/JEM-X/ISGRI spectrum could be well fit with a cut- 
off power-law model (F = 0.5 + 0.2, E Mi = (5.7 + 0.6) keV, 
£cut = (10.8 + 0.6) keV), b ut some significant residuals were 
left above 10 keV. Following Tiirl er et al. we showed in 

Sect l2.2l that the fit could be improved by adding a cyclotron ab- 
sorption feature with a centroid energy at (22 +1) keV. The fea- 
ture was not significantly detected during the latest stages of the 
outburst (from MJD 55 991 onward), most likely owing to the 
relatively low S/N of the data. As the centroid energy of the line 
lies in the energy range in which the JEM-X and ISGRI instru- 
ment responses overlap, we caution that at present it is impos- 
sible to firmly exclude that some calibration uncertainties might 
have affected this detection. Future observations are thus needed 
to establish its presence. This feature would however not be un- 
expected in the case of an X-ray pulsar. 

A hard spectrum with T < 0.6 is sometimes observe d 
from sources that display cyclotron lines (Cob urn et al.l l2002). 
and could be related to the effects of cyclotron scattering 
(Ferr igno et alj 120091) . Cyclotron resonant-scattering features 
(hereafter CRSFs) might also appear in the spectra of these ob- 
jects in the form of absorption lines. The latter are caused by 
resonant scattering of photons off the electrons in Landau levels 
in the strong magnetic field of the neutron star (10 u -10 13 G). 
The centroid energy of the fundamental absorption feature is 
related to the strength of the m agnetic field in the scattering 
region by the equation (see e.g., iWasserman & Shapirol ll983) 
£ cycl ^11.6/(l+z)(B/10 12 G) keV, where B is the magnetic field 
strength in gauss and z+l=1.31 is the gravitational redshift (we 
used the canonical neutron-star mass and radius of 1 .4 M0 and 
10 km, respectively). If the cyclotron feature is confirmed in the 
future, we will be able to infer for IGR J18 179-1621 a magnetic 
field strength of ~2xl0 12 G. 



4 



E. Bozzo et al.: IGR J18179-1621 : an obscured X-ray pulsar discovered by INTEGRAL 



A similar magnetic field strength, associated with pulsations 
at ~ 12 s, would make IGR J 1 8 1 79- 1 6 21 a likely B e X-ray binary 
candidate (hereafter BeXBs; see e.g. lReigll201lL for a r eview). 
Acco rding to this interpretation, the Corbet diagram dCorbetl 
Il986h suggests for this source an orbital period in the range 
20-50 days. Despite this apparently straightforward association, 
such a relatively short orbital period makes the interpretation of 
the event recorded by INTEGRAL and Swift in terms of the 
so-ca lled BeXB "type-I" outbursts challenging (see e.g., iReigl 
1201 ll) . As the latter occur almost regularly at each periastron pas- 
sage in these systems, previous detections of IGR J18 179-1621 
should probably be expected during the long-term monitoring of 
the sky around the source performed so far with INTEGRAL. 
By using the online tool HEAVENS^, we checked that the total 
effective exposure-time at the coordinates of IGR J18 179-1621 
is 2.4 Ms for ISGRI and 410 ks for JEM-X (considering all pub- 
licly available data from 2003 February 28 to 2010 March 1). 
A search in these archival data did not result in any previous 
significant detection of the source, and we derived a 5cr upper 
limit to the source flux of 0.5 mCrab in the 17-40 keV energy 
band and 1.7 mCrab in the 3-10 keV energy band. The inter- 
pretation of the event re corded from IGR J 18179-1 62 1 as a rare 
BeXB type-II outburst (Stel la et al.lll986l) . might also be ques- 
tionable. The X-ray luminosity reached during these events is 
typically on the order of Lx-10 38 erg s~ l , and would thus imply 
an unrealistically large distance for IGR J 1 8 1 79- 1 62 1 (» 1 kpc, 
using the peak flux measured in Sect. 12.21 ). If the event reported 
here is the first episode of intense X-ray activity ever displayed 
by IGR J18 179-1621, this source could be a further example of 
those "dormant" X-ray binaries that spent most of their time life 
in a quiescent state before suddenly changing int o bursting X-ray 
sources (see e.g. the ca se of RX J0440.9 +443 1 : lTsvgankov et alj 
l2^TllUsuiet"atll20T2h . 

At odds with typical BeXBs, IGR J18179-1621 appears to 
be characterized by a peculiar high absorption column-density. 
A value of jVh-10 23 cirT 2 is more commonly observed in the 
so-calle d highly obscure d high-mass X-ray binaries (HMXBs; 
see e.g. lChatv et al.l[2010L for a review). A case of a highly ab- 
sorbed HMXB displaying a cycl otron scattering fea ture is the 
persistent source GX 301-2 (see ISuchv et al.ll2012l and refer- 
ences therein). In these sources, the neutron star is accreting 
matter from the intense wind of its massive companion giving 
rise to a fairly persistent X-ray luminosity in the range 10 34 - 
10 36 erg s (the exact value depends mainly on the strength of 
the wind and on the orbital separation between the supergiant 
and the neutron star). Bright short flares, lasting a few hours and 
reaching a peak flux ~ 10-50 higher than the persistent emission, 
are sporadically observed from these objects and ascribed to 
episodes of accretion from high density "clumps" of stellar wind 
material din't Zandl2005HBozzo et aDl2Q08t iRrev kenboh m et al] 
120081; iNegueruelal 120101 : IBozzo et al.1l201 ll) . The dynamic range 
of the X-ray flux measured from IGR J18 179-1621 (taking into 
account the upper limit derived above), might thus still be 
marginally compatible with that measured for the highly ob- 
scured HMXBs, although the scale time of the outburst decay 
inferred from the Swift monitoring (~22 days) seems hardly rec- 
oncilable with that expected for an episode of clumpy wind ac- 
cretion. 

Even though some of the results from the timing analysis 
of the XRT data have to be interpreted with caution (see lAl . in 
Sect. [2] we reported the possible detection of a spin-period varia- 
tion from the source. If this is interpreted in terms of an accretion 
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torque, the average spin-up and spin-down rates inferred from a 
linear fit to the measurements before and after MJD 55 989.13 
(see Fig. [3} would be on the order of ~10~ 7 s/s, i.e. compatible 
with thos e measured in a number of accreting high -mass X-ray 
binaries (Bild sten et al.|[l997t iFerrigno et al.l 20071 see also the 
case of GX 1-4). We note that similar changes in the spin period 
of the source might also be due to the orbital modulation, but the 
limited number of currently available estimates of P sp i„ do not al- 
low us to investigate further this possibility. Future observations 
of the source in outburst are needed to strengthen and confirm 
any detection of spin-period variations in IGR J18179-1621. 

As an alternative interpretation, we also suggest that 
IGR J18 179-1621 might be another example of those peculiar 
binaries displaying a complex X-ray variability with interme- 
diate properties in-between those of young wind-accreting sys- 
tems and more evolved disk accretors. The latter can display pro- 
longed periods of quiescence and undergo moderate to bright 
outbursts ow ing to instabilities in the accretion disk (see e.g., 
Lasota 2001, and references therein). Among the known pecu- 
liar disk-accreting systems, a particularly relevant case here is 
Her X-l, which is also known t o display a cyclot ron absorp- 
tion line at ~40 keV (see e.g., IVasco et aD 1201 ll and refer- 
ences therein). Spectroscopic observations of IGR J18179-1621 
in other energy domains (e.g. IR) might help in clarifying the 
nature of its companion star and unveil the real nature of the 
mechanism regulating the X-ray variability displayed by this 
new INTEGRAL transient. 
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Appendix A: Spin-period determination 

The uncertainty and the significance of the pulse period derived 
using the Z 2 statistics in each of the XRT observations reported 
in Sect. 12.21 could be reliably estimated only by performing 
Monte Carlo simulations, as the individual determinations ob- 
tained by scanning frequencies from a single event file are not 
independent. Furthermore, we verified that the formal signifi- 
cance of the Z 2 peak is not a good estimate of the reliability of 
the measurement for the Swift observations. The latter comprise 
several relatively short (%\ ks) snapshots containing only a few 
hundreds photons each and, as a result, the corresponding peri- 
odogram has the shape of a wide-peak (the diffraction figure of 
the window function) with evident interference fringes produced 
by the fragmentation of the observations (see Fig. lA.U . 

To perform the Monte Carlo simulations, we determined in 
each observation the best pulse period and approximated the cor- 
responding pulse profile by using its first two Fourier compo- 
nents (a check was performed a posteriori to verify that this de- 
scription is accurate in each case by inspecting the single pulse 
profiles; see Fig. [2j». Event files with an average event rate equal 
to that measured from the source were then simulated in the good 
time intervals (GTIs) of the observations by introducing two si- 
nusoidal modulations at the estimated spin period and its first 
harmonic. The amplitudes of the modulations were assumed to 
be the same as those measured from the Fourier decomposition 
of the real source pulse profiles. In each of the simulated event 
files, the pulse period was then measured with the Z 2 -statistics 
method and recorded. We report the results of this analysis in 
Fig.EH 

The spin periods determined from the simulated files cluster 
around the real value measured in the corresponding observa- 
tions and follow roughly a Gaussian distribution. The shape of 
the Gaussian could in principle be used to estimate in each case 
the uncertainty in the spin-period determination, but this proce- 
dure is complicated in all cases by a certain number of "outliers". 
The number of these points, whose "wrong" value of the spin 
period is due to the sporadic displacement of the peak in the Z 2 - 
statistics periodogram at the moment of the best period determi- 
nation, is larger for shorter observations and obser vations char- 
acteriz ed by a smaller number of source events. In lBozzo et alJ 
(1201 lh . we argued that the number of outliers with respect to 
the total number of measurements could give an estimate of the 
reliability (significance) of the spin-period detection in an obser- 
vation. 

To identify the outliers in the present case, we: (i) collected 
all the spin periods determined from the 90% of the simula- 
tions that gave periods closer to those measured from the real 
data ("central spin periods"), (ii) estimated the average spin pe- 
riod and the variance cr 2 of the corresponding Gaussian distribu- 
tion, and (iii) identified all the remaining realizations that gave a 
spin period differing by more that 2.6<x from the averaged one. 
We rejected as "not-reliable" the XRT observations in which the 
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Fig.A.l. Periodograms obtained from the Z 2 statistics applied 
to the XRT observations reported in Sect. I2.2I The OBSIDs are 
indicated by the two last digits in each panel. 

number of outliers were found to significantly exceed (99% c.L, 
taking into account also the intrinsic uncertaint y in cr and the 
Poissonian nature of the key variables; see e.g., Sheskrn1 l2007l) 
the value expected from the Gaussian distribution of the simu- 
lated central spin periods. 

As a final remark, we note that, in all cases, the timing anal- 
ysis was performed on XRT event files that were uncorrected 
for pile-up (see sect. [2}, as the latter is known to affect mostly 
the spectral energy distribution of photons recorded from the 
soured We checked a posteriori that an analysis similar to 
that described above performed on the pile-up corrected XRT 
event files of the reliable spin-period determinations in observa- 
tions 00032293001 and 00032293002 would give values fully 
in agreement (to within the uncertainties) with those reported in 
Table For the observations 00032293004 and 00032293008, 
the lower flux and/or number of counts led, in contrast to an un- 
reliable spin-period determination once the corresponding data 
are corrected for pile-up. As the effect of pile-up on these obser- 
vations cannot be checked further owing to the limited statistics, 
the corresponding values of the spin periods reported in Table |4] 
should be interpreted with caution. 



4 See http://www.swift.ac.uk/analysis/xrt/pileup.php 
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Fig. A.2. Spin periods derived from the simulated event files for 
each of the XRT observation reported in Sect. 12.21 The OBSIDs 
are indicated with the two last digits for each panel. 
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